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Physical unitarity for a massive Yang-Mills theory without the Higgs field: 

(I) A perturbative treatment 

Kei-Ichi KondoQ Kenta Suzuki, Hitoshi Fukamachi, Shogo Nishino, and Toru Shinohara 
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In a series of papers, we examine the physical unitarity in a massive Yang-Mills theory without 
the Higgs field in which the color gauge symmetry is not spontaneously broken and kept intact. For 
this purpose, we use a new framework proposed in the previous paper ^lj based on a nonperturbative 
construction of a non-Abelian field describing a massive spin-one vector boson field, which enables 
us to perform the perturbative and non-perturbative studies on the physical unitarity. In this 
paper, we present a new perturbative treatment for the physical unitarity after giving the general 
properties of the massive Yang-Mills theory. Then we reproduce the violation of physical unitarity 
in a transparent way. This paper is a preliminary work to the subsequent papers in which we present 
a nonperturbative framework to propose a possible scenario of restoring the physical unitarity in 
the Curci-Ferrari model. 
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I. INTRODUCTION 

In the Yang-Mills theory |2J and quantum chromo- 
dynamics (QCD) for strong interactions, both the renor- 
malizability and the physical unitarity are satisfied, as 
demonstrated first in Q. Moreover, it is also known that 
the massive Yang-Mills theory in which the local gauge 
invariance is spontaneously broken by the Higgs field [4] 
and the gauge field acquires the mass through the Higgs 
mechanism by absorbing the Nambu-Goldstone particle 
associated with the spontaneous symmetry breakdown 
satisfies both the renormalizability and the physical uni- 
tarity In other words, both the renormalizability and 
the physical unitarity survive the spontaneous breaking 
of the gauge symmetry. 

It is a long-standing problem [rjl-[l6j to clarify whether 
it is possible or not to construct a massive Yang-Mills 
model blessed with both the physical unitarity and the 
renormalizability without the Higgs fields, in which the 
local gauge symmetry is not spontaneously broken. Here 
the Lagrangian is assumed to be written in polynomi- 
als of the fields (we exclude the nonpolynomial type [l3[ 
from our discussions). Such a model is anxious to un- 
derstand the mass gap and confinement caused by the 
strong interactions jl7H2lj , since there does not exist the 
Higgs field and the color gauge symmetry is kept intact 
in QCD. Indeed, there are continued attempts to look 
for an alternative way to describe massive non-Abelian 
gauge fields without the Higgs field ill. However, all 
these efforts were unsuccessful in coping with both renor- 
malizability and unitarity very well: In all the models 
proposed so far for the massive Yang-Mills theory with- 
out the Higgs fields, it seems that the renormalizability 
and the physical unitarity are not compatible with each 
other, although there are some models which satisfy ei- 
ther the renormalizability or the physical unitarity. See 
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[iH for reviews and [l6| for later developments. 

For this purpose, we start once again from the Curci- 
Ferrari (CF) model @, which is a massive extension of 
the massless Yang-Mills theory in the most general renor- 
malizable gauge having both the BRST and anti-BRST 
symmetries [22|. In preceding studies for the CF model 
|9l-[l2|. the CF model is proved to be renormalizable, 
whereas the CF model has been concluded to violate 
physical unitary [Iol - [T2| . However, the preceding stud- 
ies are restricted to considerations in the perturbation 
theory. We need a nonperturbative framework to draw a 
definite conclusion to this issue. 

In a previous paper therefore, we have presented 
a nonperturbative construction of a massive Yang-Mills 
field which describes a non-Abelian massive spin-one 
vector boson with the correct physical degrees of freedom 
without the Higgs field |4J . This is achieved by finding a 
nonlinear but local transformation from the original fields 
in the CF model to the physical massive vector field 
which is invariant under the modified BRST and anti- 
BRST transformation. As an application, we have writ- 
ten down a local mass term for the Yang-Mills field and 
a dimension- two condensate, which are exactly invari- 
ant under the modified BRST transformation, Lorentz 
transformation and color rotation. The resulting mas- 
sive Yang-Mills model is regarded as a low-energy effec- 
tive theory of QCD, which enables us to understand the 
decoupling solution fl8j characterizing the deep infrared 
regime responsible for color confinement (23l . |24| . 

In a series of papers, we give the perturbative and non- 
perturbative studies on the physical unitarity [2oT | in the 
massive Yang-Mills theory constructed in the previous 
paper [l|. In the ordinary massless Yang- Mills theory, 
the physical unitarity is a first step of understanding color 
confinement [23j : In the intermediate state, the contri- 
butions from the unphysical gauge modes, i.e., the longi- 
tudinal and scalar modes are exactly canceled by those 
of the ghost and antighost, which is a special case of the 
quartet mechanism [26}. We clarify how the situation 
changes in the massive case. Moreover, we clarify the 
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reason for failures of the preceding attempts from our 
viewpoint. 

This paper is the first of the planned papers for dis- 
cussing the perturbative and non-perturbative physical 
unitarity in the massive Yang-Mills theory without the 
Higgs field. In this paper we present a new perturba- 
tive treatment for the physical unitarity after reviewing 
the general properties of the massive Yang-Mills theory. 
Then we reproduce the violation of physical unitarity in 
a transparent way. In subsequent papers, we present a 
nonperturbative framework to discuss a possible scenario 
of restoring the physical unitarity in the massive Yang- 
Mills theory. 

Finally, we mention the difference between the unitar- 
ity and physical unitarity of the scattering matrix from 
our point of view. For the tree-level scattering amplitude 
between two longitudinally polarized vector bosons, it is 
known (27l - [29j that the scattering probability as a func- 
tion of the energy E becomes greater than 1 above a crit- 
ical value E c , since the amplitude grows with the energy 
E like g 2 E 2 /M 2 where M is the mass of the vector bo- 
son and g is the coupling constant for the self-interactions 
among vector bosons. This implies that the perturbative 
unitarity breaks down in high-energy E > E c . There- 
fore, for the perturbative unitarity to be satisfied, the 
energy must be restricted to low-energy E < E c , which 
is called the unitarity bound. (The violation of the uni- 
tarity condition for the scattering amplitude in high en- 
ergy is understood from the Nambu-Goldstone equiva- 
lence theorem [29j | and the low-energy theorem.) In the 
Higgs sector of weak interactions in the standard model, 
the Higgs particle exists and the exchange of the Higgs 
particle affects the amplitude so that the amplitude ap- 
proaches a constant at energies far above the Higgs pole. 
Consequently, the Higgs mass must be less than an up- 
per bound. If such new physical degrees of freedom do 
not exist, this behavior is not modified and the unitarity 
violation in high energy cannot be avoided in the mas- 
sive Yang-Mills theory, since the Nakanishi-Lautrup aux- 
iliary field can be integrated out and the ghosts can play 
no role in the tree level amplitude. In our works, we re- 
gard the CF model as a low-energy effective theory of the 
Yang-Mills theory to be valid in the region E < E c for 
discussing color confinement. We restrict our examina- 
tion on the physical unitarity to sufficiently low energy 
region below a few GeV to evade the unitarity violation 
and consider only the physical unitarity, i.e., unphysical 
mode cancellation in our papers. Therefore, the well- 
known fact about the unitarity violation in the above 
does not contradict with our research on the physical 
unitarity. 

This paper is organized as follows. 

In section II, we introduce a massive Yang-Mills the- 
ory without the Higgs field and define the CF model as 
a special case. The CF model is not invariant under the 
usual BRST and anti-BRST transformations. However, 
the CF model can be made invariant by modifying the 
BRST and anti-BRST transformations. The cost of in- 



troducing the mass term is the violation of nilpotency of 
the modified BRST and anti-BRST transformations. We 
point out an important fact that even the modified BRST 
(anti-BRST) invariant quantity depends on a parameter 
P in the M case. This should be compared with 
the M = case, in which (3 is a gauge-fixing parameter 
and the BRST-invariant quantity does not depend on (3, 
which means that the physics does not depend on j3 in 
the M = case. This is not the case for M ^ 0. 

In section III, we summarize the result obtained in the 
previous paper [l| on a nonperturbative construction of 
a non-Abelian massive Yang-Mills field J^, under the re- 
quirements which guarantee (i) the modified BRST (and 
anti-BRST) invariance, (ii) correct degrees of freedom 
for describing a massive spin-one particle, and (iii) the 
expected transformation rule under color rotation. We 
write down the massive vector field explicitly in terms of 
the original Yang-Mills field, the Faddeev-Popov (FP) 
ghost field, antighost field and the Nakanishi-Lautrup 
(NL) field in the CF model. 

In section IV, we give a perturbative framework of the 
CF model in terms of new field variable J^. We give the 
Feynman rules up to the order g. 

In section V, we check the physical unitarity in the 
massless Yang-Mills theory. Using a simple example, 
it is demonstrated in the lowest order of perturbation 
theory that the physical unitarity follows from the can- 
cellation among unphysical modes: the longitudinal and 
scalar modes of the Yang-Mills field together with the FP 
ghost and antighost. 

In section VI, we review a conventional argument for 
the violation of physical unitarity in the massive Yang- 
Mills theory without the Higgs field. Using a simple ex- 
ample corresponding to the previous section, we show 
that the violation of physical unitarity follows from the 
incomplete cancellation among unphysical modes: the 
scalar mode with the FP ghost and antighost. 

In section VII, we begin with a new analysis on the 
physical unitarity of the CF model based on a novel 
framework using the field given in section III. In this 
section, we give a new perturbative analysis using the 
result of section V. We confirm that the physical unitar- 
ity is indeed violated in the CF model in the framework 
of the perturbation theory in the coupling constant. It 
is easily seen that the violation of physical unitarity fol- 
lows from the incomplete cancellation among unphysical 
modes: the NL field (corresponding to the scalar mode) 
with the FP ghost and antighost. We discuss how to 
avoid the violation of physical unitarity within the per- 
turbative framework. 

In the final section, we summarize the results and men- 
tion the perspective on the next work. 

In Appendix A, we calculate the Jacobian associated 
with the change of variables from the original CF model 
to the new theory written in terms of new variables. In 
Appendix B, the Feynman rules are given up to the next 
order g 2 , with which we supplement the results of section 
V. 
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II. THE CURCI-FERRARI MODEL AND THE 
MODIFIED BRST TRANSFORMATION 

In order to look for a candidate of the massive Yang- 
Mills theory without the Higgs field, we start from the 
usual massless Yang-Mills theory in the most general 
Lorentz gauge formulated in a manifestly Lorentz covari- 
ant way. The total Lagrangian density is written in terms 
of the Yang-Mills field st^ the Faddeev-Popov (FP) 
ghost field the antighost field and the Nakanishi- 
Lautrup (NL) field JV . As a candidate of the massive 
Yang-Mills theory without the Higgs field, we add the 
"mass term" J<? m : 



"*mYM 



1 



se> 



GF+FP 



-^GF+FP 



(la) 
(lb) 



a 
2 2 

+ JV ■ &*sfp - ^gJT ■ (iff x tf) 

+ itf-d> 1 

+ ^(,yf ? -,yF + ,yV-,yV) + ^,yy-,yV, (lc) 
1 



M 2 ^ ■ ^ + 0M 2 i c € ■ 



(Id) 



where a and are parameters corresponding to 
the gauge- fixing parameters in the M —> limit, 

9y\^\€{x) := d^{x) + gsi{x) x <tf(x), and 



JV := -JY + gitf x , 



(2) 



The a = case is the Curci-Ferrari (CF) model with 
the coupling constant g, the mass parameter M and the 
parameter 0. In the Abelian limit with vanishing struc- 
ture constants f ABC = 0, the FP ghosts decouple and 
the CF model reduces to the Nakanishi model [30| . 

In what follows, we restrict our considerations to the 
a = case. In the a — case, Jz?ym + -£?gf+fp is con- 
structed so as to be invariant under both the usual BRST 
transformation: 



8j* m (x) = %[s/]<tf(x) 
8tf{x) = -¥S{x) x 



x) = ijV[x) 
8jY{x) = Q 

and anti-BRST transformation: 



8^(x) = (x) 
M(x) = -% c iix) x 'rf(x) 
8^{x) = ijV(x) 
8,jV{x) = {) 



(3) 



(4) 



Indeed, it is checked that 

8^ YM = 0, <SjS?gf+fp = 0, (5) 
S^ym = 0, <5_£?gf+fp = 0. (6) 

This is not the case for the mass term Jz? m , i.e., 

sse m o. (7) 

Even in the presence of the mass term ^f m , however, 
the total Lagrangian -2^ym can be made invariant by 
modifying the BRST transformation Q: ^b R st = 
with a Grassmannian number A and 



'S'^(x) = %[*/}tf(x) 
S'tf(x) = -§tf(x) x ¥(x) 
6"?ix) = ijV(x) 
8'Jf{x) = M 2c €(x) 



(8) 



The modified BRST transformation deforms the BRST 
transformation of the NL field and reduces to the usual 
BRST transformation in the limit M — > 0. It should be 
remarked that S'^f^ YM = follows from 



while 



o — <5 (-S?gf+fp + se m ), 



S'j?m 0, <5'J2?GF+FP ¥= 0. 



(9) 



(10) 



Similarly, the total action is invariant under a modified 
anti-BRST transformation 5' defined by 



5"r?(x) = -%V(x) x ¥(x) 
8'^(x) = iJ?(x) 
8',jV{x) = -M 2c £(x) 



(11) 



which reduces to the usual anti-BRST transformation in 
the limit M — > 0. It is sometimes useful to give another 
form: 



8'jV{x) =gjT{x) x <f( x ) - MV(i), 
8'jV(x) =gJV{x) x <if(x) + M 2c i{x). 



(12) 



Moreover, the path-integral integration measure 
V^V c €V c iVjV is invariant under the modified BRST 
transformation. Indeed, it has been shown in [l[ that the 
Jacobian associated to the change of integration variables 
$(x) -» $'(x) = $(x) + \8'$(x) for the integration mea- 
sure is equal to one. 

However, the modified BRST transformation violates 
the nilpotency when M^0: 



8'8'^(x) = 0, 
S'8"€{x) = 0, 

8'8 u i(x) = i8',yV(x) = iM 2 tf{x) ^ 0, 

8'8'jV{x) = M 2 8'tf{x) = -M 2 ^(x) x V{x) ^ 0. 



(13) 
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The nilpotency is violated also for the modified anti- 
BRST transformation when M ^ 0: 



8'8'^{x) = 0, 
8'8"rf(x) = 0, 
8'8"i?(x) = i8'J?{x) 

8'5'j?{x) = -M 2 8 lc i{x) = M 2 ^(x) x <§{x) ^ 0. 



-iM 2 tf(x) ^ 0, 
'.gc 

2 ' 



(14) 



In the limit M — > 0, the modified BRST and anti- 
BRST transformations reduce to the usual BRST and 
anti-BRST transformations and become nilpotent. 

Moreover, it is checked that the modified BRST 
and modified anti-BRST transformations no longer anti- 
commute in the M/0 case: 



{8\8'}^(x) = 0, 
{S',S'}f(x) = -iM 2r ^(x), 
{8',8'}Y(x) = iM 2c i{x), 
{{8>,8'}JS(x) = 0. 



(15) 



In the limit M — > 0, the anticommutativity is recovered, 
{<5',5'}^0. 

Let W be the generating functional of the connected 
Green functions defined from the vacuum functional Z[J] 
with the source J for an operator & as a functional of $: 



iW[J] 



(16) 

Then the derivative of W with respect to B is given by 



dw[j] = i din z[j] = i x dz\j\ = / asK™ 

dp i 83 i 1 1 d/3 \ d/3 



Z[J}:= j WBOTP/ 

x exp \ iS^y M + i / d D x.J(x) ■ &(x) 



where 



(17) 



mYM 



d/3 



= d 



(18) 

This follows from the fact that JzfcF+FP + -S?m is written 
as Q 

+ !lM 2 i c i-<tf+ ^M 2 ^ -j**. (19) 

If we require the modified BRST invariance of the vac- 
uum: 



QbrstIO) = 0, 
we find the (3 dependence of W[J]: 



dW[J] 

d/3 



d D x-M 2 {itf(x) 



(20) 



(x))j^0. (21) 



Therefore, for M ^ 0, W[J] depends on the parameter 
8, This result should be compared with the M = case, 
in which B is a gauge fixing parameter and hence W[J) 
should not depend on 8. In the M = case, any choice of 
8 gives the same W[J]. However, this is not the case for 
M 7^ 0. The 8 dependence of the CF model was pointed 
out also in [33f using different arguments. 



III. DEFINING A MASSIVE YANG-MILLS 
FIELD 

We require the following properties to construct a non- 
Abelian massive spin-one vector boson field J^(x) in a 
non-perturbative way: 

(i) has the modified BRST-invariance (off mass 
shell): 



8'j^ = 0. 

(ii) is divergenceless (on mass shell): 



d^JXTf, = 0. 



(22) 



(23) 



(iii) obeys the adjoint transformation under the 
color rotation: 

,^{x) -> UJtT^U- 1 , U = exp[ie A Q A ], (24) 



which has the infinitesimal version: 

<5J^ t (x) = £ x JX^(x). 



(25) 



The field Jt^ is identified with the non-Abelian version of 
the physical massive vector field with spin one, as assured 
by the above properties. Here (i) guarantees that be- 
long to the physical field creating a physical state with 
positive norm, (ii) guarantees that Jt^ have the correct 
degrees of freedom as a massive spin-one particle, i.e., 
three in the four-dimensional spacetime, i.e., two trans- 
verse and one longitudinal modes, excluding one scalar 
mode, (iii) guarantees that J^, obey the same transfor- 
mation rule as that of the original gauge field srf^ 

We observe that the total Lagrangian of the CF model 
is invariant under the (infinitesimal) global gauge 
transformation or color rotation defined by 



(26) 



5<f>(x) := [e c iQ c , $(ar)] = e x $(x), 

for $ = ^ lll ^,^, c ^, 
8tp(x) := [e c iQ c , tp(x)] = -ietp(x), (27) 



where tp is a matter field, which is also written as 
6<f> A ( x ) = f ABC e B $ c {x), 

Sip a (x) = -is A (T A ) a b tp b (x) = -ie A (T A ip) a , (28) 

where the conserved Noether charge Q A := J d 3 x ^ikx 
obtained from the color current ^^loi * s cane d the color 
charge and is equal to the generator of the color rotation. 
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It has been shown [l| that such a field Jif^ is obtained 
by a nonlinear but local transformation from the original 
fields si^ , and =yf of the CF model: 

- M -2 ^^ - gM~ 2 ^ x 

+ gM~ 2 d^ x i^ + g 2 M~ 2 (^ x x i c #. 

(29) 

In the Abelian limit or the lowest order of the coupling 
constant g, Jt^ reduces to the Proca field for massive 
vector: 

^ sty. - -^pd^ ■= U„. (30) 

It should be remarked that Un is invariant under the 
Abelian version of the modified BRST, but it is not in- 
variant under the non- Abelian modified BRST transfor- 
mation. 

The new field Jt^ is converted to a simple form: 

= + ( 31 ) 

It has been explicitly shown in [l| that the field Jt^ de- 
fined by (|29[) or (|3~Tj) satisfies all the above properties. 
The field plays the role of the non- Abelian massive 
vector field and is identified with a non- Abelian version 
of the spin-one massive vector field. The equation (|2"9")l 
gives a transformation from srf^ JV 1 ^ and c <§ to J^. 

As an immediate application of the above result, we 
can construct a mass term which is invariant simultane- 
ously under the modified BRST transformation, Lorentz 
transformation and color rotation: 

^M 2 J^(x)- Jtr»(x). (32) 

This can be useful as a regularization scheme for avoid- 
ing infrared divergences in non-Abelian gauge theories. 
Moreover, we can obtain a dimension-two condensate 
which is modified BRST invariant, Lorentz-invariant and 
color-singlet: 

(J^(z) • JT»{x)). (33) 

This dimension-two condensate is off-shell (modified) 
BRST invariant and should be compared with the 
dimension- two condensate proposed in [3~ll |32| : 

(^(x)-^(x)+l3tf(x)-tf(x)y (34) 

which is only on-shell BRST invariant. 

The original CF Lagrangian -S^m^,^, ^,J^} is 
written in terms of stf^,^,^ and . The new theory 
is specified by i^Ym[J^„ "W, ^Y] written in terms of 
J^, % ^ and J/ with the symmetry: 

(5'.^{x) = 

j 5'tf{x) = UT{x) K ' 

\&.jV{x)=M 2 ^[x) 



IV. PERTURB ATIVE FRAMEWORK OF THE 
MASSIVE YANG-MILLS THEORY 

The eq. (|29|) gives a transformation from stf^,^,^ 
and JV to J^. In order to write explicitly the non- 
Abelian massive Yang-Mills theory without the Higgs 
field, we rewrite the original Lagrangian written in terms 
of srffj,, JY , ^ and ^ into the new Lagrangian written in 
term of Jfc^jY, ^ and . For this purpose, we need 
the inverse transformation of (f2U)) . namely srf^ as a func- 
tion of J^, , and c <a . But the inverse transformation 
cannot be given in a closed form, since (|29|) is a nonlin- 
ear transformation. In order to perform the perturbative 
calculation, it is sufficient to know the order by order ex- 
pression of the inverse transformation. By using (|29p in 
the form: 

dp =<%n + M~ 2 dyjV + gM' 2 ^ x JY 

- gM~ 2 d^ x tfi- g 2 M- 2 {srf fi x tf) x iff, (36) 

where the right-hand side contains the order g°,g l and 
g 2 terms. By iterative procedures, we obtain 

^ =J^ + M~ 2 d^ - gM- 2 itf x 

+ gM- 2 ^ x,yV + gM-^d^jV x,yV + 0(g 2 ). 

(37) 

The propagator of Jf^ is obtained from the order g° 
terms, i.e., s/^ = Jf^ + j^d^^f, as in the Abelian case, 
and hence the propagator is not modified from the Proca 
case. However, the vertex functions among J^, "iff, and 
,yV are modified in the massive theory are modified by 
the order g 1 terms from those of srf^ ^, ^ and JV in the 
original theory. This fact has been overlooked in all the 
preceding studies. The preceding works [6j-|l2j are based 
on the observation that the vertices in the massive case 
are the same as the massless case. However, this ob- 
servation is correct only if the relationship between the 
original field srf^ to the massive field J^, is linear as in 
the Abelian case. This is not the case in the non-Abelian 
case, as our construction of the massive field clearly 
shows. The vertices in the massive theory is modified in 
terms of the vector field J^, in addition to the FP ghost, 
antighost and the NL field. For the correct identification 
of the massive vector field J^j, one needs the FP ghost, 
antighost and the NL field, in addition to the original 
Yang-Mills field s^^. This could be a loophole of avoid- 
ing the results of the preceding analyses for the violation 
of physical unitarity. It is checked if this is true or not. 

To check the physical unitarity in the nontrivial lowest 
order in the perturbation with respect to the coupling 
constant g, the original Lagrangian Jz?^y M is rewritten 
order by order of the coupling constant g as follows. 

if CF = if +^ 1 + O(.g 2 ), (38) 
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k 

to ^ wvwwv je? 

k 

(b) y I > <gB 

k 

(c) jyA .jyB 

FIG. 1: Feynman rules for propagators: (a) massive vector 
propagator, (b) FP ghost propagator, (c) NL field propagator. 

j% = - \{o^ - d v x p f + ^m 2 x p • x» 

- • fl M Sf + f3M 2 i^ ■ (39) 
" jjf*^ ' x 

+ ^( d ^ - ■ x V* 

- id vc i x fl"Sf ) 

- g X p • (iSf x S'V) + 5 JT P • (i<9 M<r # x Sf ) 

+ ' x ^) - fff ^ ' («f x ). 

(40) 

The Feynman rules are given up to three-point vertices 
of 0(g) as follows. 
Propagators (Fig. Q} 

(a) massive vector propagator (X A (k) X B (— k)) 

-W^TTe(^-^) sAB < 41 » 

(b) FP ghost propagator ( c tf A (k) c i B \-k)) 

(c) NL (auxiliary) field propagator (X A (k)X s (— k)) 

iM 2 - R , N 

S AB (43) 

k 2 -(3M 2 + ie [ ' 

Three-point vertices (Fig. [2|) 

(d) XXX vertex (X A {k x )X B (k 2 )X p c '(fc 3 )) 

-igf ABC V^ p (ki,k 2 ,k 3 ), (44) 



(<l) 




wvww ^ 



(f) 



9f s • 



c#A 



hi 



(g) 



FIG. 2: Feynman rules for 3-point vertex functions: (d) 
XXX vertex, (e) X^ vertex, (f) XXX vertex, (g) 
vertex. 



with 

V fivp (h,k2,fo) :=(h - k 2 ) p g l _ iv + (fc 2 - k 3 )^gup 

+ (*3 - fci)^5pM- (45) 

(e) JT^*? vertex (^(h)^ 13 {k 2 )X p c (k 3 )) 

- i gf 4sc {A^ 2 [fc lp (fc 2 ■ fc 3 ) - fc 2p (fc! • fc 3 )] 

+ ki p - k 2p }. (46) 

(f) XXX vertex (^ A (fci)^ s (fc 2 )JT p c (fc 3 )) 

i5 Af- 2 / ABC {Af- 2 [fc lp (fc 2 • fc 3 ) - k 2p (k x ■ fc 3 )] 

+ kip - fc 2p }. (47) 

(g) Xtf^ vertex ( < ^ A (fci) c # ,B (fc 2 )^ /C (fc 3 )) 

-*3-/3/2]. (48) 

These rules should be compared with the Feynman 
rules for the original Lagrangian: 
Propagators (Fig. [3]) 

(a) vector propagator (g/ A (k)&/ B (— k)) 
i 



k 2 — M 2 + is 



9»v 
kirk. 



kXv \ cAB 



k 2 + is 



k 2 - f3M 2 + is k 2 + ie 



V 6 AB = D AB (k). (49) 
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k 

( a ) K A / vwvwv\y K B 

k 

(b) ■/ I ► y; B 

k 

(c) ^4 ^/B 

k 

( d ) 'WW ^ s 

FIG. 3: Feynman rules for propagators: (a) massive vector 
propagator, (b) FP ghost propagator, (c) NL field propagator, 
(d) cross propagator. 

w ;:*\AA/w £/ p c 

cgB 

FIG. 4: Feynman rules for 3-point vertex functions: (e) 
sisrfsrf vertex, (f) s4 c g < i vertex. 

(b) FP ghost propagator {^ A {k) c i B {-k)) 

(c) NL field propagator {jV A {k)^ B \-k)) 

(d) the cross propagator {s/ A {k).jV B {— k)) 

k 2 -/3M 2 +ie ■ V ; 

Three-point vertices (Fig. 0]) 

(e) ,06^ vertex (^(fciXffej^fe)) 

-^/ ASC ^p(fci,fc 2 ,fc 3 ), (53) 

(f) vertex ( < # A (fci)^ B (fc 2 )^' p c '(fc3)) 

igfABchp- (54) 



From the relation, 



(55) 



we find 

(^(aj)^(y)) 

= (PW + M- 2 d^{x)){J(r v {y) + M~ 2 d^{y))) 

= {,^(x),je v {y)) + M- 2 {,^{x)d u ^{y)) 

+ M- 2 {d^jr{x)Jf v {y)) + M- A {d^{x)d v jV{y)), 

(56) 

which has the Fourier transform: 

=(jr/(fc)jr/(-fc)> + ik u M- 2 {jr A {-k)^ B {k)) 
-ik ll M- 2 {j^ A {k)jr v B {-k)) 



+ M- 4 k f ,k u (,yK A (k),yY B (-k)) 



--id 



AB 



M 2 



M 2 



k 2 - M 2 + ie 



j3M 2 



(57) 



where we have used (|41|>. (|49p and the fact that there are 
no mixing propagators: 

(je A {-k)^ B {k)) = o = (,A- A (k\je B {-k)). (58) 

Thus the original gluon 2-point function or propagator is 
decomposed into the spin-1 and spin-0 parts. 

Another expression for the propagator is the manifestly 
(power-counting) renormalizable form: 



«(fc)< B (-fc)) 
9nv 



--iS 



AB 



k 2 +ie 







W+T. 



k 2 - M 2 + ie k 2 - /3M 2 
For higher orders, see Appendix B. 



(59) 



V. PHYSICAL UNITARITY OF MASSLESS 
YANG-MILLS THEORY 

The S'-matrix or the scattering operator S is unitary: 

1 = 5+5 = 55 f . (60) 

This means that for any (initial) state £ V and any 
(final) state vj/ a £ V, the following relation holds: 



(tf 6 |* ) =<* 6 |5+5|tf a ) = £ (* 6 |5t|$ n )($ n |5|* a ), 

(61) 



which is obtained by inserting the complete set of states 
{$„} in the total state space V: 1 = J2<t> ev \ ®n)(®n\- 

On the other hand, the physical unitarity of the 5- 
matrix means that the S matrix is unitary on the physical 
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subspace V p h ys : for any state physical state *b <= 

Vphys i 

(* 6 |*.> = J2 (*b\S*\$ n )(* n \S\*a), (62) 
*„GV phys 

where the physical subspace V p h ys is defined by 

Vphys := {|phys) e V; (phys|phys) > 0} C V. (63) 

The unitarity of the S-matrix is rewritten in terms of 
the scattering amplitude defined by 



into the relation: 



S = l + iT, 



i(T-T f ) = TTl 



(64) 



(65) 



Then the unitarity relation reads: for any state ty ai £ 
V, 

Im(* b |T|* Q ) :=l«* 6 |T|* a ) - (* 6 |Tt|* a )) 

=\ E <*&m*„><*»|:rt|* a >. (66) 

*„ev 

On the other hand, the physical unitarity requires that 
for any physical state ^ a ,^b & Vphys, only the physical 
states contribute to the intermediate states: 





FIG. 5: In the massless case M = 0, the diagrams contribut- 
ing to the amplitude T(.e/ — ► s/) to the order g 2 are given 
by (a) vector boson loop, (b) ghost-antighost loop, (c) boson 
tadpole. 




FIG. 6: In the massless case M = 0, the physical unitarity 
of the amplitude T(T ± — > T ± ) to the order g 2 is checked 
according to the Cutkosky rule using the diagrams: (a) vector 
boson loop, (b) ghost-antighost loop. 



The imaginary part of the diagram Fig. H][a) with a 
iduon loop is given by 



Im<* fe |T|* a } = - ]T (* 6 |T|$ n )<<& n |Tt|* a ). (67) 

#„£ Vphys 

In other words, the physical unitarity in gauge theories 
states that all the unphysical modes cancel in the inter- 
mediate states. The imaginary part is calculated by the 
Cutkosky cutting rule [34j . 

As a simple example, we consider a one-particle scat- 
tering, i.e., a scattering process si — > in which the 
initial state is a massless gluon and the final state is also 
a massless gluon in the massless Yang-Mills theory as the 
M = case of the CF model. The Feynman rules in the 
M = case have been given in [32j]. In the lowest order 
of the perturbation theory, the Feynman graphs of this 
process are given by Fig. [5j Each diagram has one closed 
loop. The initial state and the final state of a gluon are 
specified by the polarization vectors: 



(k) 



.(o-)a' 



(k) 



(68) 



For this process up to the order g 2 , we wish to check the 
physical unitarity relation for the scattering amplitude 
T. By applying the Cutkosky rule [34| to Fig. [5J we 
find that the imaginary part of the scattering amplitude 
T(T ± — > T^) from a transverse mode to a transverse 
mode is given by Fig. [5] Here it should be remarked 
that the tadpole diagram does not have the imaginary 
part. 



{ig)fABcVvvct {h ,k 2 ,k 3 ) 



1 

x {-ig)fABC'V IJ ,' t/ > a '(-ki, -fc 2 , -fc 3 ) 

x 2ne(kl)8(kl)g^'2^e{kl)8(kl) 9 ^'e a {k 3 )e a ' (-fc 3 ), 

(69) 



where we have adopted the Feynman gauge j3 = 1 for the 
gluon propagator and the factor 1/2 is the symmetrical 
factor due to two identical particles. This is written as 



Ig^fABcfABC'T^gW <T T M ^^(fc?)<5(fc?)0(fc 2 °) ( 5(fc 2 2 ), 



where we have defined 



(70) 



:= V PlUa {k 1 ,k 2 ,k 3 )e a (k 3 ), 
Z>„ :=^w(-Aa,-fc2,-fc3)£ Q '(-fc3). (71) 



In the massless case M = 0, the physical unitarity 
requires the imaginary part of the sum of the second- 
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order diagrams to be equal to 

Ig^fABcfABCiT^P^'P^'T^v, 

x 6(kl)5{k\)e{kl)5{kl) 
= ^(ig)fABcV f _ ll y a (k 1 ,k 2 ,k 3 ) 

x (-ig)fABC'V fl >„> a >(-ki, -k 2 , -fc 3 ) 

x 27r0(fc?)$(fc?)P w * ' 2Ti9{k£,)5(kl)P uv ' ' e a (k 3 )e a ' (-k 3 ) 

1 2 

= 2 E ( l .9/^cV«(fci,fc2,fc 3 ))^(fc 1 )4(fc 2 )e Q (fc 3 ) 

a,6=l 

x (igfABcVuw (fci , fc 2 , fc 3 ))*er' (fciK"' (k)e ' 

x 2 7 r6'(fcJ)5(fc 2 )27r6'(^)5(^), (72) 

which is obtained from (|70[) by the replacement: 

2 

a=l 
2 

s -%p-'=^4(fe) £ *"'(fc 2 ), (73) 

b=l 

where P correspond to the two transverse polarization 
states for the massless spin-one modes T^. 
By using the decomposition: 

ff""' = -i 3 ""' + Q vv ' > Q vv ' = (kfkZ + W)/(*a ■ h), 

(74) 

the difference between = (170)) and (|72")l is calculated 
from 

gW' g vv ' — piw' p vv ' 

= - pw' Q vv ' - P vv ' Q^' + Q^' Q uv ', (75) 

where Q are rewritten using the polarization vectors for 
the longitudinal (L) and the scalar (S) modes: 

QT 1 = e£ (h )ef (hi) + eg (fci )e*"' (A*), 

Q""' = el(k 2 )ef(k 2 )+sUk2)ef(k 2 ). (76) 

By using the relationships: 

T^{k x )el{k 2 )=T» v el{k 1 )el{ki) 
I (-i)k la e a (k 3 ) (a = S) 







(<r = T,L) ' 



(77) 



and 



= Uk la ,e a '(-k 3 ) ((7 = 5) 
10 (<r = T,L)' 




T • 

+ 1-1)— Q 



FIG. 7: In the massless case M = 0, mode cancellations 
occur to ensure the physical unitarity for the one-particle 
amplitude T(T ± — > T ) for the transverse mode T to the 
order <? 2 . In the amplitude, two diagrams (a) from the lon- 
gitudinal mode L and the scalar mode S are cancelled by a 
ghost-antighost c € diagram (b). 



we find that the first and second term in the right-hand 
side of (|75p give vanishing contributions. 

The relationship ([77)1 is derived as follows. First, we 
find 

V^aiki, k 2 , k 3 )e u L (k 2 ) 
=V^ a (ki,k 2 ,k 3 )k 2 

= {ki - k 2 ) a k 2 ^ + (k 2 - k 3 )^k 2a + (k 3 ■ k 2 - ki ■ k 2 )g a ^ 
= - hahfi + k 3a k 3fll (79) 

where we have used k\ — and k 3 = for the massless 
on-shell momenta. Then we have 

iWfci, fc 2 , fc 3 )^(fciK(fe)e Q (fc3) 
= - k^ik^k^ih) + k 3a e a {k 3 )k 3 ^ a {ki), (80) 

where k 3a e a (k 3 ) = for the physical transverse mode, 
while fci M £y(fci) = 0, ki^e^{k\) = —iki^k^ = for mass- 
less on-shell k\, and kx^e^ikx) = i for massless on-shell 
h. 

Thus, only the last term in (|75[) gives a non- vanishing 
contribution: 



4:g 2 Tr 2 fABcfABC>ki a e a (k 3 )k la ,E a '(-k 3 ) 
x 0(fc?)<5(fc 2 Mfc 2 o )<5(fc 2 2 ). (81) 



(78) 



This difference is exactly provided by the imaginary 
part of the second diagram of Fig^b) with a ghost loop: 

(-l){-igfABcki a )(igfABC>ki a >)£ a (k 3 )e a (-k 3 ) 

x 2Tr9(k ( l)5(kl)2ir9(k 2 > )6(kl) 

= - 4:g 2 ir 2 f AB cfABC'ki a e a (k 3 )k la ,£ a ' \-k 3 ) 

x 0(*J)$(*?)0(*§)<y(*§), (82) 

where we have used the original gluon-ghost-antighost 
vertex (given in [32j at £ = 1/2 corresponding to a = 0): 
\igfABc{ki a -k 2a ) = igfABc{ki a + ^k 3a ) and the prop- 
erty of polarization vectors: k 3a £ a (k 3 ) = to obtain the 
first expression. Thus, the unitarity relation is satisfied 
in the massless case M = 0. 

The physical unitarity is assured by mode cancella- 
tions. See Fig. [7] The contributions from the longitu- 
dinal mode L and the scalar mode S are cancelled by a 
ghost-antighost ^ one. 
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VI. PHYSICAL NON-UNITARITY OF MASSIVE 
YANG-MILLS THEORY 

In this section, we reproduce the violation of physical 
unitarity in the massive Yang-Mills theory without the 
Higgs field based on the conventional argument, which 
shows the utility of the massive field J(f in discussing the 
physical unitarity of the CF model in the next section. 

In order to see a difference between massless gauge 
theory (M — 0) and massive vector theory [M ^ 0), 
we consider the one-particle scattering, i.e., a scattering 
process U — > U in which the initial state is a massive 
vector boson and the final state is also a massive vector 
boson. Here is defined by 



1 1 



M 2 



M 2 



(83) 



In the lowest order of the perturbation theory, the 
Feynman diagrams of this process are given by the same 
graphs as those in Fig. [5] where the propagators are re- 
placed by the massive ones and the vertex functions are 
unchanged. Therefore, the imaginary part is given by the 
diagrams of Fig. [8] 

The imaginary part of the diagram Fig. HJa) with a 
loop of massive vector boson is 

\ {-ig) fABcVfj.ua (h , h, k 3 ) 

x (ig)fABC'V^ 1/ ' a '(-k 1 ,-k 2 ,-k s ) 
x 2TT0(k° 1 )6{k 2 1 - M 2 )g^' 

x 2^9{k° 2 )5{k 2 ~ M 2 )g^'s a (k 3 )s a '(-k 3 ), (84) 

where the factor 1 /2 is the symmetrical factor due to two 
identical particles. The imaginary part (|84|) is written as 



Ig 2 '* 2 fABcfABC'Tfj V g^ g vv 



x Q{k\)8{k\ - 
where we have defined 



M 2 )6(k%)5(kl - M 



(85) 



Tf W ■■= Vf ll , a (ki,k 2 ,k 3 )e a (k 3 ), 



T„ 



V^w(-*i,-*2,-*3)e a (-*&)• (86) 



In the massive case the physical unitarity re- 

quires the imaginary part of the second-order diagram to 
be equal to 



1g 2 K 2 f A BcfABc>T^P^ P vv Tfj, v , 
x e{kl)5{kl - M 2 )6(k° 2 )5(kl - M 2 ) 

= \ {igfABcTfjvE^ (ki )e^{k 2 )) 

x iigfABC-T^^' (fci)e?'(fca))* 

x 2ir0(k° 1 )6(k 2 1 - M 2 )2-K9{kl)5(k 2 - M 2 ), (87) 



M f 




FIG. 8: In the massive case M 7^ 0, the physical unitarity 
of the amplitude T(U U) for the physical spin-one vector 
mode U to the order g 2 is checked according to the Cutkosky 
rule using the diagrams: (a) vector boson loop, (b) ghost- 
antighost loop. The physical unitarity is violated due to the 
incomplete cancellation. 




5 ~ Q 



FIG. 9: In the massive case M 7^ 0, the incomplete mode 
cancellation among unphysical modes to the order g 2 prevents 
us from ensuring the physical unitarity for the one-particle 
amplitude T(U — > U) for the physical spin-one vector mode 
U. In the amplitude, a diagram (a) from the scalar mode S is 
overcancelled by a ghost-antighost diagram (b), leaving one 
half of (b) nonvanishing. 



which is obtained from (|85[) by the replacement: 

g^' ^ P^' ^e^ih^ih), 

9 VU ' p uv ' '=£^)(M4>)- 



Using the decomposition for M 7^ 0: 



M 2 



I"'' +^p-, (89) 



where P selects three polarization states (T + ,T~ , L) for 
the massive spin-one boson U, the difference between 
and (|57|) is calculated from 



(90) 



j K 1 "'2 puv' K l K l 1 K l K l ft 2 K 2 



pfj,fj. 



M 2 



M 2 M 2 M 2 



The contribution from the first term of the right-hand 
side of to {S5) is given by 



2g 2 j^TT 2 f AB cfABC'Tfj^(k 1 )k^ (fci)A^T^/ 
x e{k\[)5(k 2 - M 2 )e{k° 2 )5{k 2 2 - M 2 ). (91) 
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This is zero, since 
T^{kx)k v 2 

^e';(k 1 )V fwa (k 1 ,k 2 ,k 3 )k^e a (k 3 ) 

= effa^-kipkia + /s 3 M fc 3a + kfg^a - klg fia ]s a (k 3 ) 

= -h^ik^k^ih) + k 3 , t e$(k 1 )k 3a e a (k 3 ) 



0, 



(92) 



where we have used k\ = k 2 and k^s^Jk) = 0. Sim- 
ilarly, the contribution from the second term of (f90|) is 
vanishing. The third term of (|90|) gives 



k^ky k? k v 

Using the property: 

Vfcf k v 2 = -k\k la e a {k 3 ) + kx ■ k 3 k 3a e a (k 3 ) 

= ~M 2 k la e a (k 3 ), (94) 

we obtain the difference: 

2g 2 TT 2 fABcfABC>ki a e a {k 3 )k la ,e a '{~k 3 ) 
x 9{kl)5{kl - M 2 )6{k° 2 )5(k 2 - M 2 ). (95) 

This difference must be provided by the imaginary part 
of the diagram Fig. [5Jb) with a loop of massive ghost 
which is given by 

(-l)e a {k 3 )e a ' (~k 3 ){-igf A Bcki a )(igfABC>(-k 2a >)) 
x 2ir0(k ( l)5(kl - M 2 )27T0(k^)S(kj - M 2 ) 
= - 4:g 2 TT 2 f AB cfABC'ki a e a {k 3 )k la ,e a '(-k 3 ) 

x 6(k%)6(k 2 - M 2 )9{kl)5{k 2 - M 2 ), (96) 

where we have used k 2 = — k\ — k 3 and k 3a e a {k 3 ) = 0. 
The ghost contribution (|96|) is precisely of the same form 
as (|95|) and comes with the opposite sign. However, the 
massive vector contribution (p?5"j) cancels against half the 
massive ghost contribution ([96)). See Figj9] 

Thus, it is found that there is a discrete difference be- 
tween massless theories and massive theories. This means 
that the massless theory cannot be obtained as a limiting 
case of the massive theory. The origin of the difference 
goes back to the difference between the sum over polar- 
izations. In the next section, we re-examine the physical 
unitarity in the the CF model based on our method. 



VII. PERTURB ATIVE ANALYSIS OF 
PHYSICAL UNITARITY IN THE CF MODEL 

In order to re-examine the physical unitarity in the 
massive Yang-Mills theory, i.e., the CF model, we con- 
sider the simplest case of the one-particle amplitude 
T(J(f — > J(f) in the perturbation theory, as considered 




X ; «■ «■■••. x xi ,Jr ■ yr \x 



FIG. 10: In the massive case the physical unitarity 

of the amplitude T(,3tf — > ,X) to the order g 2 is checked ac- 
cording to the Cutkosky rule using the diagrams: (a) massive 
vector boson, (b) ghost-antighost, (c) NL field. 



According to the Cutkosky rules, the physical unitarity 
up to the order g 2 is checked by calculating the diagrams 
in Fig. [TU] with one closed loop. 

The imaginary part for the diagram Fig. HOf a) is 

(a)S£ := Q) (ig)fAB C V^ a (k u k 2 ,k 3 ) 

x {-ig)fA>Bic>V»'»' a '(-k u -k 2 ,-k 3 ) 
x 2-K0{kl)8 AA ' ' I^{k{)2-Ke{kl)5 BB ' I vv ,{k 2 



where we have defined 



I» v {k) :=5{k 2 - M 2 ) g^ 



~M 2 



(97) 



(98) 



Remember that the three polarization vectors effl (k)(j — 
1,2,3) for the spin-one massive vector field obey the 
relation: 



3 

E 



^\k 2 )e^*{k 2 ) = -g vul 



M 2 ' 

k 2v k 2v > 
M 2 ■ 



(99) 



Therefore, (a) is a physical contribution coming from the 
spin 1 massive vector boson J(f . The physical unitarity 
requires that the contributions other than (a), i.e., the 
contributions (b) and (c) from the unphysical fields c £ ', & 
and ,JV are cancelled in the same order of the coupling. 
Therefore, we consider the contributions from unphysical 
fields. 

The imaginary part for the diagram Fig. UOf b) with a 
closed loop of the massive FP ghost and antighost reads 

(b) cc : 

■■={-l)(-ig)fABc{M~ 2 V a (ki,k 2 , k 3 ) + k la - k 2a } 
x igfA'B'c{M~ 2 V a '(ki, fc 2 , k 3 ) + kia> - k 2a '} 
x 2it6{k° 1 )5{kl - PM 2 )2i:6{k^)5{kl - j3M 2 ) 

= - 45 2 7r 2 Jabc f abc {M~ 2 V a {ki,k 2 ,k 3 ) + ki a - k 2a } 
{a -> a'}6(k^)6{kl - f3M 2 )e{k° 2 )5{kl - pM 2 ). (100) 
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M 4 



j + (-ll~-H^) 




FIG. 11: In the massive case M^O, incomplete mode can- 
cellations violates the physical unitarity for the one-particle 
amplitude to the order g 2 . 



The imaginary part of the diagram (c) with a closed 
loop of the NL field reads 



(c) 



cc 

aa' 



\ J ^'^M2' fABC'{M 2 V a (k 1 ,k 2 ,k 3 ) + ki a - k 2a } 
x ij^fA>Bic>{M- 2 V a ,(k 1 ,k 2 ,k 3 ) + k la , - k 2a ,} 



L0\xn„2 _ p M 2\ x AA' A/T 2 

2n6{kl)5{kl - (3M 2 )5 BB ' 
=2g 2 TT 2 f ABC fABC'{M~ 2 V a {k ll k 2 , fc 3 ) + ki a - k 2a } 
{a -> a'}8{kl)5(k 2 - pM 2 )B{k° 2 )5{k 2 - j3M 2 ). (101) 

Hence, adding the NL loop (c) to the ghost loop (b) yields 
the half of (b): (6)+(c) = \{b), since (c) = There- 
fore, we have shown the CF model does not satisfy the 
physical unitarity for M/0 case independently from /3 
in the perturbation theory. The incomplete cancellations 
of the unphysical modes against the physical unitarity is 
summarized in Fig. 111! This result should be compared 
with the massless case where the physical modes in the 
massless case M = are the two transverse parts. 

In the massless case M = 0, for the amplitude T — > T 
of the transverse mode T , two diagrams (a) from the 
longitudinal mode L and the scalar mode S are cancelled 
by a ghost-antighost f ' diagram (b). The NL mode 
,vY is non-propagating in the massless case and does not 
contribute to this cancellation. 

In the massive case M ^ 0, for the amplitude J(f — > Jif 
of a physical particle JfT, a diagram (a) from the scalar 
mode S (— the NL mode JY) is not sufficient to can- 
cel the ghost-antighost 'if, diagram (b) . An additional 
bosonic contribution is necessary to realize the complete 
cancellation. 

In the massless case M = 0, the physical modes are 
given by two transverse modes T^. Then the two un- 
physical modes in the gluon, i.e., the longitudinal mode 
L and the scalar mode S are cancelled by a ghost and 
antighost 'if, 'if. Two bosonic modes are exactly cancelled 
by two fermionic (anti-commuting) modes. It should be 



remarked that the scalar mode is identified with the NL 
mode JY , which is non-propagating in the massless case. 



T+ T~ 



L, S(= JY : non - propagating); <i , (102) 



In the massive case M ^ 0, the physical modes are 
given by a longitudinal and two transverse modes. A re- 
maining unphysical mode, i.e., a scalar mode is not suf- 
ficient to cancel the ghost and antighost contributions. 
Therefore, the elementary fields in the original action of 
the CF model are not sufficient to respect the physical 
unitarity. An additional (bosonic) mode must be sup- 
plied by some reasons. (Note that the NL field jV is 
propagating in the massive case and therefore the NL 
mode is expected to play the important role in the can- 
cellation in the massive case, in sharp contrast to the 
massless case. However, the NL mode is identical to the 
scalar mode on mass shell and hence cannot be counted 
as another independent field.) 



T+, T~,L ,S = JV; 



(103) 



Finally, we discuss how to avoid the violation of phys- 
ical unitarity. The violation of the physical unitarity is 
avoided by restricting the relevant energy to the low- 
energy region such that the ghost and antighost pair can- 
not be created. This can be done by adjusting the pa- 
rameter in the CF model. Since the ghost and antighost 
has y/]3M, the allowed region is 

E < 2V/3M, E 2 < 4/3M 2 . (104) 

A shortcoming of this scenario is that j3 — is not al- 
lowed to maintain physical unitarity, since the results of 
numerical simulations on the lattice are available only in 
this case ,8 = 0. 

At first glance, the cancellation of unphysical modes 
works well even in the massive case by using the argument 
similar to that done in the gauge-Higgs model with the 
renormalizable gauge. The pole masses of unphysical 
fields ffjff, JY are the same: 



mlg 



f3M 



j2 



(105) 



In the limit (3 — > oo, unphysical fields f, f, JY decou- 
ple from the theory, leaving the physical field Jff^ in the 
theory. Therefore, it seems that the physical unitarity 
holds even in the massive case. However, this is not the 
case, as we have shown in the above. What is wrong in 
this argument? This argument is based on the fact that 
P is a gauge fixing parameter and the physics does not 
depend on this parameter /3, which is indeed shown in 
the massless case M = 0. However, even the BRST in- 
variant quantities depend on /3 for M ^ 0. Therefore, 
the physics depends on /3 for M ^ 0, and the physics for 
P — > oo is different from that for /3 < oo. In this sense, 
the above result, i.e., violation of physical unitarity in 
the massive case does not contradict with this argument. 

In order to maintain the physical unitarity in the mas- 
sive Yang-Mills theory without the Higgs field, therefore, 
we need a non-perturbative approach, which will be given 
in subsequent papers in preparation. 
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VIII. CONCLUSION AND DISCUSSION 

In order to understand color confinement in QCD in 
the light of recent developments, we have considered a 
"massive" Yang-Mills model without the Higgs field, es- 
pecially, the CF model, since the CF model is regarded 
as a good low-energy effective theory of QCD and it is 
much simpler than the refined Gribov-Zwanziger model, 
see e.g., [l|. 

We have examined the physical unitarity of the CF 
model which is known to be renormahzable. For this 
purpose, we have used the field Jtf^{x) with the following 
properties: 

(i) Jff is invariant under an extended BRST transfor- 
mation, S'J(f fl (x) = (off shell). 

(ii) Jff is divergenceless, d' J 'Jt fl (x) = (on shell). 

(iii) transforms according to the adjoint representa- 
tion under color rotation. J^(x) — > U 'J^(x)L r_1 

(iv) is invariant under the FP conjugation. 

J^(x) -> J^(x) 

Thus, we have identified ,j4f A with a physical and massive 
vector field with correct degrees of freedom as a non- 
Abelian spin-one massive vector boson. is obtained 
by a nonlinear but local transformation from the original 
fields in the CF model. 

We have checked in a new perturbative treatment 
whether or not the CF model satisfies the physical unitar- 
ity. Then we have confirmed the violation of the physical 
unitarity in the perturbative treatment and we have clar- 
ified the reason in the massive Yang-Mills theory without 
the Higgs field. The perturbative analysis for the physi- 
cal unitarity impose a restriction on the valid energy to- 
gether with the parameter of the CF model: E 2 < Af3M 2 
in order to confine unphysical modes (ghost, antighost, 
scalar mode). However, ft = is not allowed in this 
scenario. 

The conclusion obtained in this paper still leaves a pos- 
sibility that the nonperturbative approach can modify 
the conclusion. In a subsequent paper, indeed, we will 
propose a scenario in which the physical unitarity can be 
recovered in the CF model thanks to the FP conjuga- 
tion invariance. Indeed, we will show that the norm can- 
cellation is automatically guaranteed from the Slavnov- 
Taylor identities if the ghost-antighost bound state ex- 
ists. In this way, the physical unitarity can be recovered 
in a non-perturbative way. To show the existence of the 
bound state of ghost and antighost, the Nambu-Bethe- 
Salpeter equation is to be solved. This is a hard work to 
be tackled in subsequent papers. 

Acknowledgements: The authors would like to thank 
Professor T. Kugo for very enlightening discussion about 
some issues. This work is supported by Grant-in- Aid for 
Scientific Research (C) 24540252 from Japan Society for 
the Promotion of Science (JSPS). 



Appendix A: Change of variables 

The original theory is given by 

z mYM = J v^vtfvtfv^re iS £™i^^,v,^] (A1) 

The exact change of variables: srf^ V, JV — > 
J£^, If, ^, jV could be performed through the relation- 
ship stf — stf\J<f ill ' : € 1 c io 1 Jf\ according to 



J mYM 



apr) 



= VJfV^V^V^e lS "™^^^^\ (A2) 



where 

Cm <£, Jf\ := & K, Jf\ - i In J, 

(A3) 

where J is the Jacobian associated with the change of 
variables from stf^^ ^ ,jV in the original Lagrangian 
to J£^, "if , "if, JV in the new theory. 

We proceed to check the modified BRST (mBRST) in- 
variance of the new theory. The action is mBRST invari- 
ant by construction. We have already shown that the in- 
tegration measure 'Dsi'D'^'D'^'Djy is mBRST invariant. 
The mBRST invariance of the measure VJfV^V^VjV 
is also checked in the same way. Therefore, the Jacobian 
J must be mBRST invariant too, i.e., 5' J = 0. 

However, we do not know the exact expression of 
£/ = Jl^,^,^,^} as a functional of J^, %f, 
although we know the exact expression of — 
J^l^/fi , ff, c to, JV\ as a functional of ^ , "if, ^, jY as given 
in (|29[) . We know just the order by order relation for 
si = stf\W il , c € ^ as given in (|3"T)) . Hence we can 
calculate the Jacobian J order by order of the coupling 
constant g. Thus, the mBRST invariance of J can be 
checked order by order in the coupling constant g, al- 
though the full mBRST invariance cannot be checked 
because we do not know the exact form of J. Here it 
should be remarked that S ' Jf = M 2c € and 8 lc € = iJf 
do not change the order, while 8 lc £ = ~\gi^ x ^) and 
= ^^Is/]^ increase the order of g by one. 

In fact, the Jacobian J is calculated as follows. The 
integration measure is transformed as 



The Jacobian J is calculated as 



(A4) 



r 



=Det 
=Det 



6< B (y) 
8»JS AB 



gM -2 f ABC^C 



+ g'M~ 2 f ACE f BFE t^ u ^)5 1J (x - y) 



(A5) 
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where we have used 



Appendix B: Lagrangian and Feynman rules in the 
order g 2 



M 



±_ f AEC f EGF^F^C_ 



(A6) 



Applying the formula: 

Det(l + X) 
= exp[trln(l + X)] 



■ exp 



tr(X) - ±tr(X 2 



+ ltr(X 3 ) - itr(X 4 ) + 



(A7) 



to 



X AB = -gM- z f Aau ^ u + g z M-~ z f AU *f at *i^<tf 



-2 (ABC A/-C 



2 rACE f BFE -c^Cc/pF 



(A8) 



we hnd that the order g contribution (from the NL field) 
vanishes due to f AAC = 0: 



J- 1 =]Jcxp 

X 

1 



-g 2 M~ 4 N c ,yV ■ ,jV + g 2 M- 2 N c i't? ■ c g 



+ ^ g Hr & f ABE j BCF f CAG ^v E ,^ F ^v G 

- g 3 M- 4 f ABE f BCF f GAG ,yV E itf F tf G + 0(.g 4 ) 



(A9) 



where we have used jAecjAef = NcS cf Therefore] the 
correction from the measure to the Lagrangian density 
begins with the order g 2 . In other words, J = 1 up to 
the order g. 

This is also checked by using the relationship up to the 
order g 1(57]): 



The relation s/ = s/[J^, ft, tf, jV] up to the order g 2 
is obtained as 

2 2 



„2 



M 2 



The transformation (IB1I) leads to 



Ss/ A (x) 



(Bl) 



-_ 5 »( 5 AB +gM -2 f ABC^C 



5Jtr*(y) -" v 

-g 2 M- i f AGE f BFE ,A' c ,A' F 

- g 2 M' 2 f ACE f BFE i<i c tf F )5 D {x - y) 

+ 0(g 3 ). (B2) 



Indeed, this gives the same Jacobian as (|A9j) up to the 
order g 2 . In order to check the consistency of the Ja- 
cobian (|A9|) . therefore, we must obtain the Lagrangian 
density Jzf2 in the order g 2 . 

We must take into account the correction coming from 
the Jacobian J in obtaining the Lagrangian density ,S?2 = 
jS?2p^, ^\ m the order g 2 . Then the Lagrangian 
density _£?2 in the order g 2 is obtained as follows. 



J% = 



£ if 2 W -ilnJ( 9 2 ), 



(b) 



(a) 



„2 



2M 4 
2M 4 

-,2 



— (J^j X J^) , 



X J/) ■ (STX* X ,_yf) 



- c? 2 M- 4 (J^ x ^f) • (^Jf" x 9^) 

+ <? 2 a/- 4 (j^ x ,yf) • [ff'jf x 



g^g = S;(S AB +gM- 2 f ABG .yr G )S D ( X -y) + 0(g 2 ). 

(A10) 

Thus, the Lagrangian density = Jz?i[J£T, ^, up 
to the order g does not change and remains in the form 
(14TJ1) . and the calculation obtained using _£fi (|4U)) is not 
affected. 



j^ c ) = + fM^id"^ x ,yK) • (<9^ x d u Jf), 



2M 2 



(Jf x ,yK) 2 



2 2 

■ 9 {d^xd v ^f+ 9 



4M 8 



2Af 6 



x ,yVf, 
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(e) 



if. 



(/) 



- g 2 M~ 2 {d il d il X v x fif) • (X x if) 
+ g 2 M~ 2 (d v d ll X li x fif) • (X x Sf) 
+ g 2 M- 2 (d"J(r v x X) ' x <W 

- g 2 M- 2 (d u Jt fl x X) ' (*^ x 9u^) 
+ g 2 M- 2 (W» x JT) ■ (id^ x 9,%?) 
+ 5 2 (X x *?) • (Jf" x 

- g 2 M- i (d ti d> i Jf v x flf) • x Sf) 
+ g 2 M- i (d fl d v Jt tl x fif) • (d v JV x Sf) 
+ g 2 M~ 2 {d IJ ,J/ x %f) • (JT^ x «f) 

+ g 2 M~ 2 {^ x ) • (ifl^^f x if), 

2 

+ g 2 M- 4 {d„^Y x Jf) ■ (id^ x Sf) 

- g 2 M~ i (d fJ ,JK x .jV) • (flf x 0"Sf), 



(ft) 



+ 



2M 4 
2M 4 

„2 



(id^g x 9^ 



2M 2 



(id^ x 9^) • (i9 v <«? x d^) 



M 2 



{i<€ x d^) ■ x *f) 



4 v ; 



(B3) 



The following are the Feynman rules for the vertex 
functions of the order g 2 obtained from (|B3j) . 

(a) X X X X vertex: 



(ki)-K B {k 2 )X p c {k 3 )X a D (fc 4 )) 

= - Q 2 W ABCD 



(B4) 



with 
W, 



ABCD 



{fM j,»jj _ fA u,u» )gimgpa 
+ (f AB - CD f AD ' BC )9» p9v „ 

+ (f AC ' DB f AB ' CD )g^ P , 

jAB,CD _jABE jCDE 



(B5) 



(b) XX J/ J/ vertex: 



(X A (*i)^ a (*2)^°'(*sK D W) 



M 2 

t 

M 4 



AC.BD _|_ jAD,BC 



rAC,BD 



JJi9^ [( fc l ' fc 3 + &2 ' ^4 + fei • fc 2 )/ y 

+ (fc a • fc 4 + fc 2 • fc 3 + fci • fc 2 )/ AABC ] 



M 4 



AC.BD 



+ (hvkin + k Xv k^ + k 3u k2 fl )f 



AD.BCl 



g 2 



(B6) 



(c) XJfJfJf vertex: 
(X A (k 1 ),^ B (k 2 )^r c (k 3 ),yV D (k i )) 



M 6 



(B7) 



+ (k 1 -k 3 )(f AD > CB k 2 , + f AB > CD k 4 ,) 
+ {k 1 -k,)(f AC > DB k 2 , + f AB - DC k^) 

(d) JfJfjYjy vertex: 
(^(hy^ik^ih^ih)} 

a 2 

-{^{(ki-k2 + h-k i )(f AC < BD + f AD > BC ) 
+ (k 1 -k 3 + k 2 -k i )(f AB - CD + f AD ' CB ) 
+ (kt -k i + k 2 - k 3 )(f AB - DC + f AC ' DB )}. (B8) 

(e) XX<€<i vertex: 

= - M 2 (X A (k 1 )X B (k 2 )^K c (k 3 )^ D (k i )). (B9) 

(f) JT.yfW vertex: 



' M 2 



«3 M 



/ 



-r^[^-(fci-A a )^]/ AC ' M . 
(g) JfJf^ vertex: 



(BIO) 



M 4 



[(fci ■ fc 4 ) + (fc 2 • fc 3 ) - (fci • fc 3 ) - (k 2 ■ k 4 )]f 



AB,CD 



+ ^Pi ■ fc 4 )(fe • fc 3 ) - (fcr • fc 3 )(fc 2 • fc 4 )]/ 

(h) f € c i t i c i vertex: 

(^(A: 1 )^ B (fc 2 )^ c (fc 3 )^ D (fc 4 )) 

_^j-AB,CD _|_ jAD,BC\ 



AB,CD 



(Bll) 



- ^2(^2 • *0 - -^j(A 1 • /.v.! i/.-2 • A-i) 



?1 
M 4 ' 



[(fci • A; 4 + fe 2 • fc 3 )/ 



AB.CD 



+ (fci • fc 2 + fc 3 • fc 4 )/ 



+ 



A/ 4 



[(fci • fc 4 )(fc 2 • fes)/ 



ifi.CD 



+ (fcl-fc 2 )(fc 3 -fc 4 )/ 



(B12) 
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